Photon correlations in multi-mode waveguides 
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We consider the propagation of classical and non-classical light in multi-mode optical waveg- 
uides. We focus on the evolution of the few-photon correlation functions, which, much like the 
light-intensity distribution in such systems, evolve in a periodic manner, culminating in the 're- 
vival' of the initial correlation pattern at the end of each period. It is found that when the input 
state possesses non trivial symmetries, the correlation revival period can be longer than that of 
the intensity, and thus the same intensity pattern can display different correlation patterns. We 
experimentally demonstrate this effect for classical, pseudo-thermal light, and compare the results 
with the predictions for non-classical, quantum light. 

PACS numbers: 42.50.Ar, 42.79.Gn 



Periodic recurrences, or 'revivals', of an initial condi- 
tion are ageneral phenomenon common to many physical 
systems [l| . It was experimentally observed for electrons 
in Rydberg atoms Q, atomic [H, H| and molecular [E| 
matter waves, Bose-Einstein condensates ((j, and many 
more. In optics, this effect dates back to the observa- 
tions of periodic imaging of diffraction gratings by Tal- 
bot in 1836 Q, which were explained in terms of Fresnel 
diffraction by Rayleigh in 1881 [8j. The common feature 
of all these systems is that their relevant eigen-modes 
(i.e. those initially populated) accumulate phase in rates 
that are all, either exactly or approximately, integer mul- 
tiples of the same number. Thus, there exists a point at 
which the phases accumulated by the eigen-modes are 
all integer multiples of 27r, and the initial wavefunction 
is reconstructed. This effect then repeats itself period- 
ically. Furthermore, at rational fractions of the revival 
period, 'fractional revivals' occur: the initial wavefunc- 
tion reappears in multiple, shifted copies of itself. In 
applied optics, these phenomena, collectively known as 
multi-mode interference (MMI), are utilized for multiple 
beam splitting and combining [9|, [To| . Recently [ll| , such 
MMI-based multi-port integrated devices were shown to 
function as multiple beam splitters also at the single- 
photon level, displaying quantum interference effects. 

The revival of the initial intensity pattern in opti- 
cal systems through MMI is well described within wave 
optics, and do not require the quantization of the op- 
tical field. Quantum-optical effects, such as Hanbury 
Brown and Twiss [HI, or Hong-Ou-Mandel interfer- 
ence [13], affecting the relative distribution of single 
photons, can be revealed only through measurement of 
photon-correlation functions. Second-order correlation 
functions have been recently investigated for diffraction 
gratings 

dun, 

and it was shown that the recurrence 
period of the second order spatial correlation function is 
twice larger than that of the light intensity. As MMI 
systems show periodical recurrences of the light inten- 
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sity, much like diffraction gratings, and as they can be 
rather ea sily miniaturized and integrated in photonic cir- 
cuits [H, [TU, ^ is interesting to study the behavior of 
correlation functions in these systems. Two questions 
immediately arise: Do second-order correlations in MMI 
systems have a double period in respect to the total inten- 
sity, as in diffraction gratings? and how do higher-order 
correlations behave? 

Here we theoretically and experimentally explore few- 
photon correlations in MMI systems. As a simple but 
rather flexible model experimental system we use a two- 
mirror planar multi-mode waveguide. 

In this system, a general initial wavefunction recurs ev- 
ery zq = SD 2 /X lUIlcj, where D is the distance between 
the mirrors and A is the wavelength of the incident light. 
As the intensity distribution and the correlation func- 
tions of all orders are functions of the wavefunction, they 
too recur at this distance. If the input state possesses 
only trivial symmetry, meaning that its wavefunction is 
an eigenstate only of the identity transformation, zo is the 
shortest recurrence period. Thus, for such a 'completely 
asymmetric' initial wavefunction, the intensity and all 
correlations have the same recurrence period, zq. 

This situation changes when the initial wavefunction 
has a non-trivial symmetry. Such a case is that of two 
uncorrelated, identical beams entering the waveguide at 
locations symmetric about its center. This case applies to 
a quantum input state containing a single photon in each 
beam, as well as to two thermal input beams. The total 
intensity distribution along the waveguide for such initial 
states is composed of the incoherent sum of the intensity 
distributions originating from each beam separately. Due 
to the initial symmetry, the two distributions are mirror 
images of one another. Thus, for all planes where these 
distributions are composed of at most two lobes, their 
sum always results in two identical lobes. Most general! 
such planes repeat every Talbot distance, zt=zo/4: 
LLOj , as can be seen in the numerical calculation presented 
in Fig. QJa). 

The only nontrivial transformation keeping an uncor- 
related symmetric input state (either quantum or classi- 
cal) unchanged is the symmetric exchange of photons be- 
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tween the two sides of the waveguide. Such an exchange 
occurs during the evolution through the waveguide ev- 
ery z /2=2z T . Thus, the second order correlation 
function, determined by the two-photon part of the state, 
recurs every 2zt- This period is twice larger than the re- 
currence period of the corresponding total intensity. Two 
different correlations patterns can thus exist for the same 
intensity distribution. For the separated quantum input 
state, as shown in the calculation presented in Fig. QJb), 
the two photons can be either separated between the two 
beams, or bunched at the same beam. 
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FIG. 1: (a) Calculated intensity pattern for two beams prop- 
agating through a multimode waveguide of width D along the 
z direction. Light (dark) color represents low (high) intensity. 
Periodic revivals of the initial intensity pattern are seen at in- 
teger multiples of the Talbot distance zt- (b) Calculated cor- 
relation functions for two photons entering the waveguide one 
on each beam. Here the revival period is 2zt- (c) Schematics 
of the experimental system. SLM - spatial light modulator; 
MMW - multi-mode waveguide; CCD - charge coupled device. 

This principle can be extended to any case of N mi- 
correlated identical beams entering the waveguide at the 
points x n = (2n — 1 — N)D/2N. As for the case of two 
uncorrelated symmetric beams discussed above, the only 
non-trivial symmetry of the N^* 1 order correlation func- 
tion is in the symmetric exchange of photons between the 
two sides of the waveguide. The period of this correla- 
tion function is thus 2zt for all N. However, the inten- 
sity distribution recurs every time a single input beam 
divides into up to N beams, a situation occurring every 
zt/N [10]. These cases thus show a correlation recurrence 
period 2N times longer than that of the total intensity. 

Experimentally, we demonstrate these effects up to 
the 3 r< ^ order correlation function for up to 3 beams of 
thermal light. Our experimental setup is described in 
Fig. Die). Its main part is a waveguide made of two par- 
allel planar metallic mirrors, placed on translation stages 
for adjusting the distance between them - the waveguide's 
width. Two or three beams of pseudo-thermal light are 
focused on the input facet of the waveguide, and the two 
or three photon correlation function is measured at its 
output facet. This is done for different widths, D, con- 
trolling the recurrence distance, zq? f° r the fixed waveg- 
uide length, L. 

The widths used in the experiment are between 57 



and 74 /mi. The wavelength of the incident light is 
532 nm. The waveguide thus supports over one hundred 
modes. The incident light is focused on the entrance to 
the waveguide by a cylindrical lens. This is done in or- 
der to minimize the diffraction in the direction parallel 
to the mirrors. The focusing lens used was of NA^O.l, 
creating a spot about 5/im in width. Therefore, only the 
lowest ~10% of the waveguide modes are excited, and 
accurate MMI imaging is ensured [9]. Due to the finite 
conduction of the mirrors, the waveguide is slightly bire- 
fringent. To avoid any complication that may arise, the 
incident light was polarized perpendicularly to the mir- 
rors, and thus only the TM modes were excited. In order 
to create multiple spots at the entrance to the waveguide, 
the incoming light was first passed through an amplitude 
and phase spatial light modulator (SLM) used as a vari- 
able transmission grating. A second cylindrical lens was 
placed after the waveguide in order to create a magnified 
image of its output facet on a CCD camera. 

Thermal light was simulated by repeating the experi- 
ment for many different relative phases between the in- 
cident beams, controlled by the SLM. For each phase re- 
alization, the spatially resolved correlation function was 
produced from the recorded intensity, and all such single 
phase correlation functions were summed up to give the 
thermal light correlation function [l7| . 
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FIG. 2: Calculated and measured spatial intensity correla- 
tion maps for two beams at symmetric input locations. The 
first panel on the left of each row shows a schematic illustra- 
tion of the waveguide settings. Black circles denote the posi- 
tions of the input (left) and output (right) beams. The sec- 
ond and third panels present calculated correlation maps for 
quantum and thermal input states, respectively. The dotted 
squares mark the locations of the mirrors. The fourth panel 
presents the measured correlation map for pseudo-thermal in- 
put in a waveguide with L—4.85 cm, and the measured in- 
tensity distributions, for 3 waveguide settings, (a) Imaging, 
D=D Q = y/\Ljl=h7^m. (b) Equal two -way beam-splitting, 
D=66/xra. (c) Unequal two-way beam-splitting, D=12\im. 
Here the effect is specific to input beams located at zbD/4. 

Fig. [2] presents calculated and measured intensity cor- 
relation functions between two beams, introduced into 
the waveguide at symmetric locations, for three differ- 
ent waveguide settings, (a)-(c), schematically presented 
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on the left panel in each row. The second panel in each 
row presents the correlation map of light at the output 
facet of the waveguide, calculated for a quantum input 
state in which there is exactly one photon in each beam: 
ip% n = 1 11), where the left (right) number counts the 
photons on the upper (lower) beam. The third (fourth) 
panel presents the calculated (measured) output correla- 
tion map for a classical, thermal input state, where all 
photon configurations have the same probability. The 
fourth panel presents also the measured total intensity 
distribution. The two-photon part of such the thermal 
state is given by the following density matrix, written in 
the basis {|20), |11), |02)}: 

i / 1 0\ 
Pt = j 2 . (1) 
4 V 1 / 

The factor of two for the |11) state is due to the indistin- 
guishability of the photons, for which this state is reached 
from two configurations. Due to this factor and the lack 
of non-diagonal terms, thermal light shows correlations 
similar to those of the |11) quantum light, though al- 
ways with a 'background' coming from the 1 20) and 1 02) 
states pjj . 

In FigEfa) the width, D, of the waveguide is set such 
that its length, L, which is fixed at 4.85 cm, w ould 
be equal to the imaging length: D=D =y^XL/S. For 
A=532 nm, Dq= 57/im. In this case, the field at the out- 
put facet is a complete reconstruction of the field at the 
input facet. This includes the reconstruction of all corre- 
lations. For a quantum two-photon input state, the pho- 
tons at the output are also completely separated ('anti- 
bunched'), and the probability of finding two photons 
on one side ('bunching') stays zero, as seen in the cal- 
culated correlation map on the second panel, where the 
'diagonal' elements vanish. For two beams of classical, 
thermal light, the photons have the same probability to 
bunch or anti-bunch, as indeed seen in both the calcu- 
lated and measured intensity correlation maps presented 
in the third and fourth panels, respectively. 

In Fig.[2fb) the waveguide width was set to Do/ ^3/4, 
meaning that L=3zo/4=3zt- At this distance, for input 
beams at symmetric locations, the waveguide acts as an 
equal two-way beam-splitter [9l4l l| . The quantum input 
state 1 11) would therefore be transformed into the path- 
entangled state, 

^ t = -L(|20) + |02)). (2) 

The probability for anti-bunching would now be zero. 
This is seen in the calculated correlation map shown on 
the second panel. Each component of the thermal state 
pin (Eq. [1]) will undergo the beam-splitter transformation 
by itself. The resulting density matrix is: 

i / 3 1 \ 
Pfut = ^[ 2 0. (3) 
8 \ 1 3 / 



In this state it is three times more probable that the pho- 
tons would bunch than that they would anti-bunch. This 
is readily seen also in the numerical calculation presented 
on the third panel, where this ratio is the ratio between 
the sum of all diagonal lobes and that of all non-diagonal 
lobes. The bunching- anti-bunching ratio extracted in the 
same way from the measured correlation map, presented 
in the fourth panel, is 2.5±0.1. The deviation from the 
theoretical value may be attributed to background com- 
ing from parasitic effects such as mode mixing induced 
by mirror imperfections. 

Fig. [2fc) presents correlation functions for a width of 
D=D /y / tjS (L=5z T /2). For input beams at ±D/4 
from the middle of the waveguide, this is the interme- 
diate case between the cases presented in Fig. [2fa) and 
in Fig. Efb). In this case, a single photon would split 
unequally between the two beams, with amplitudes pro- 
portional to cos7r/8 and sin7r/8 [10j. The output state 
of two photons entering in the separated state is now 
an equal superposition of the separated state and the 
bunched, path-entangled state, 

^ ut = -L|ll) + i(|20) + |02)). (4) 

Therefore, there is an equal probability for bunching and 
anti bunching. This is shown in the second panel of 
Fig. [2fc). The thermal state would be transformed by 
this unequal beam-splitter into the state: 

1 ( 5 iV2 1 \ 

p'out = 6 . (5) 

\ 1 iV2 5 / 

The probability ratio between bunching and anti- 
bunching is now 5/3. This is seen in both calculation 
and measurement, the measured ratio being 1.5 ±0.1. 

It should be stressed that in all the cases presented in 
Fig.[2j the output intensity distributions are similar: they 
all include two identical lobes. The difference between 
the cases is revealed only in the correlation pattern, for 
a thermal as well as for a quantum input state. Note 
that for the special case of input beams at ±D/4, the 
correlation recurrence period is four times longer than 
that of the intensity. 

We now turn to discuss experiments with three input 
beams. We consider the case of equally spaced, identical 
beams entering the waveguide at x±=— D/3, #2=0, and 
xs=D/3. For this case, each one of the beams splits into 
three every zt/3. Although each time this three- ways 
splitting is in different weights, the total intensity of all 
three beams together is the same for all these points. The 
different weights of the single beam splitting are mani- 
fested in the correlation functions, which, as discussed 
above, recur only every 2^t, a period six times longer 
than that of the total intensity. 

For example, while for a width D=Do/ ^j2/3 
(L=8zt/3), the waveguide functions as an equal three- 
way beam splitter, for a width Dq/\/7/12 (L=7zt/3) 



it functions as an un-equal three-way beam splitter [To| . 
While the total intensity distribution for these two cases 
is the same, featuring three equal lobes, the states, and 
hence the correlations, are different. For a quantum state 
of the form 1 1 1 1 ) input into the equal three-way beam- 
splitter, the output wavefunction is, 



v. 



out 



^e i2 3 (| 300) 



|030) + |003)) 



111) 



(6) 



It carries mostly three photon bunching, along with 
some complete anti-bunching, but no two-photon bunch- 
ing [l8j |. In contrast, the unequal three-way beam-splitter 
produces the wavefunction, 



-^(e*T|210> 



^ (|300) + 2|030> + |003)) 
|201) + 1 102) 



.2tt 

■ e l 3 |012) 



(7) 



which carries mostly two-photon bunching, some three 
photon bunching, and no complete anti-bunching. Since 
we would like to have a measure for three-photon bunch- 
ing or anti-bunching, we calculate and measure the third 
order intensity correlation function. 

Theory - Quantum Theory - Classical Experiment - Classical 






50 
















E 

(b) X- 

X 























00 
X 


-50 











00 










-50 





50 


-50 50 
x 2 -x 1 (urn) 


-50 50 



Fig. [3] presents third-order correlation maps for three 
beams entering a waveguide set as either one of the two 
types of three-way beam-splitters discussed above, at the 
locations —D/3, and D/3. Fig. [3]^ a) [(b)] presents the 
case of the equal [unequal] three-way beam splitter. In 
order to conveniently present the third-order correlation 
functions, their dimensionality was reduced from three 
to two. This was done by representing them as functions 
of position differences between photons instead as of ab- 
solute photon positions, by taking the average of all ele- 
ments with the same differences. Three-photon bunching 
is represented by the middle spot. Two photon bunching 
is represented by the spots on the two main axes, and 
on the main diagonal. All other spots represent com- 
plete anti-bunching. The measurement for thermal light 
(right-most panels) is compared against the calculations 
for thermal (middle panels) and quantum (left panels) 
input states. The correspondence between experiment 
and theory is clearly seen. 



In summary, it is shown that the ratio between the re- 
currence period of the correlations between photons prop- 
agating through a multi-mode waveguide and the recur- 
rence period of their total intensity distribution depends 
on the number of input beams and on the symmetries 
of the input state. While for a totally asymmetric input 
state this ratio is 1, it can increase up to 2N for N equally 
spaced, uncorrelated, identical input beams. Cases with 
identical intensity distributions can thus display differ- 
ent correlation patterns. This is shown theoretically for 
quantum light, and both theoretically and experimen- 
tally for thermal light. Multi-mode waveguides could 
thus offer a simple yet versatile tool for the engineering 
and study of nontrivial states of light. 



FIG. 3: Spatial third order intensity correlation functions for 
three beams entering three-way beam splitters, (a) [(b)] An 
equal (unequal) beam splitter, D=70/nm [D=74fim]. The x 
and y axes present position differences between photons. 
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